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Abstract 

In this paper, we use PIC code Vorpal to do the extensive simulation about the laser 
plasma accelerator in the linear, quasilinear and nonlinear regime respectively. We design 
the lOOMeV or so laser plasma accelerator ( LPA ) via Vorpal simulation. Finally, we 
discuss the application of the designed LPA in the compact light source field. 
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1 Introduction 

Since the first proposal of laser wakefield accelerator by Tajima and Dawson[Tl [2], the laser 
plasma accelerator has gone through over thirty years[3]. In laser plasma accelerator, the lon¬ 
gitudinal Langmuir wave is excited by the driving laser which passes through the plasma. The 
perturbated electron density excites the plasma wave potential which forms the accelerating 
field as in conventional accelerator [3]. With the appearance of the quasi-monoenergetic electron 
beam from laser plasma accelerator ( LPA ), there begins the new hot advanced accelerator 
concept expedition to the laser plasma accelerator. The produced beam quality of LPA de¬ 
pends on the injection method [3 El El El El la El [E] • The present main injection methods are 
the bubble mechanism[5l EJ [7], the colliding optical pulse injection[8l |9], the controlling gradi¬ 
ent methodUni EH El ESI El ESI [E] , and the external injection met hod [T7]. After the great 
achievement of the quasi-monoenergytical electron bunch during the past ten years, it seems 
that we come to the difficult situation to go further to produce the high quality electron beam. 
We need the new breakthrough coming out of our efforts. 

2 Vorpal code 

Now, PIC code technique [I8] is widely used in the accelerator simulation field. In the past, 
we have used OOPIC[19] to do the extensive simulation about the laser plasma accelerator !^ 
[22l[23j. Vorpal [20] is a fully object-oriented code including PIC and fluid simulation, and is 
widely used in the accelerator labs. We can use Vorpal to simulate the laser plasma accelerator 
in ID, 2D, and 3D geometry. Vorpal is a good simulator for the plasma based accelerator 
simulation. My impression is that Vorpal is enough for the accelerator simulationist. According 
to our experience, OOPIC and Vorpal have become the main simulators of our simulations of the 
plasma based accelerator. In this paper, we mainly discuss the simulation of LPA via Vorpal. 
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3 Simulation of Laser Plasma Accelerator via Vorpal 


In reference [2T] , we give the simply instinctive theory of the working regimes about LPA. 
According to the laser strength, the laser plasma accelerator can be defined into three regimes: 
the linear regime, the quasilinear regime, and the nonlinear regime. In case of the linear regime 
or the linear theory, ao 1, the driving laser pulse will be guided well, and the excited plasma 
wave Wakefield has fine consine structure of a few of cycles. While, in case of the nonlinear case 
or the blow-out regime, ao ^ 1. The bubble mechanism works in this regime and can produce 
the quasi-monoenergytical electron beam. This kind of injection method work stably, but the 
Wakefield wave form is not so perfect to go further to obtain the electron beam of the more better 
quality. In case of the quasilinear regime, ao ^ 1. This is the mid way to get the high quality 
electron beam. The quasilinear regime has the partial advantages of both the linear regime and 
the nonlinear regime. The bubble mechanism, the colliding optical pulse injection, the gradient 
control injection and the external injection method can work in this regime. In the past ten 
years, the bubble mechanism has made great breakthrough in the expedition to the advanced 
accelerator concept. The relative energy spread of LPA with the bubble mechanism now stay at 
the level of few or ten percent. It may be better to return to the case of relatively small ao to 
have the plasma wave wakefield of the perfect wave form/shape. 

In order to test the incorporated code Vorpal, we do the extensive one dimension, two 
dimension and three dimension simulations via Vorpal to validate the instinctive theory of the 
three kinds of regime mentioned above. We use window platform to run the code vorpal. The 
basic laser parameters are the wavelength of 800nm, the pulse length of 1.6/im . As for the 
laser strength, we set ao = 0.1, 0.5, 0.8, 0.9,1.0,1.1,1.2, 4, 10 respectively. So we can simulate the 
linear regime, the quasilinear regime and the nonlinear regime or the bubble regime. The basic 
plasma parameter is the plasma density of 6 x 10^^m“^. The simulation results are shown in 
the Figure 1-18. Figure 1-9 show the evolutions of the electric field. While Figure 10-18 show 
the evolutions of the plasma electrons. The simulation result agree the instinctive theory result 
basically. 
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Figure 1. the longitudinal electric field ( ao = 0.1 ) 
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Figure 2. the longitudinal electric field ( ao = 0.5 ) 
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Figure 3. the longitudinal electric field ( uq = 0.8 
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Figure 4. the longitudinal electric field ( ao = 0.9 
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Figure 5. the longitudinal electric field ( uq = 1.0 ) 
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Figure 6. the longitudinal electric field ( uq = 1.1 ) 




Figure 7. the longitudinal electric field { ao — 1.2 
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Figure 8. the longitudinal electric field ( uq = 4.0 





Figure 9. the longitudinal electric field ( uq = 10.0 
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Figure 10. the longitudinal phase space ( uq = 0.1 ) 



Figure 11. the longitudinal phase space ( ao = 0.5 ) 
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Figure 12. the longitudinal phase space ( ao = 0.8 ) 
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Figure 13. the longitudinal phase space ( ao = 0.9 ) 
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Figure 14. the longitudinal phase space ( uq = 1.0 ) 
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Figure 15. the longitudinal phase space ( uq = 1.1 ) 
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Figure 16. the longitudinal phase space ( uq = 1.2 ) 
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Figure 17. the longitudinal phase space ( uq = 4.0 ) 
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Figure 18. the longitudinal phase space ( ao = 10.0 ) 


4 Design of a single cell LPA 

A lOOMeF stage laser plasma accelerator design simulation is presented. The simulation 
area is 80/im x 25.6/im, also the moving window technique is used. The used simulation PIC cell 
number is 160 x 32. The plasma density is 6 x 10^^cm“^. The laser plasma interaction works 
in the bubble mechanism. The laser strength is set to be ao = 15, the dephase distance and 
the depletion distance are estimated to be several hundreds of micrometers roughly while the 
Rayleigh length is just 20/im or so, but the laser power is far bigger than the critical laser power, 
so there exists the relativistical and ponderomotive channel optical guiding phenomenon, the 
simulated interaction range is less than the dephase and the depletion distance, the simulation 
distance is several times of the Rayleigh length. The estimated maximal energy agrees with the 
simulation result. Figure 19-21 summary the design simulation result. 


1x10 


5x10 


-5x10 


-1X10 



.OOOE+007 

.eooE+oos 

.40GE + GGa 

.200E + 008 

.OOOE+OOS 

.aoGE+GGa 

.G6GE + GG3 




9 




7x10"^ 4xin"'' Gx10"^ Rx10 


^-5 


-5 


\-5 


7 



















Figure 19. the plasma density 



Figure 20. the longitudinal electricfield ( uq = 15 ) 
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Figure 21. the logitudinal phase space ( uq = 15 ) 


5 Design for Thomson backscattering light source 

The electron beam drived x-ray light source can be defined into two classes: coherent light 
source and incoherent light source. In case of high energy electron beam, free electron laser is 
the coherent light source, while the synchrotron radiation is the incoherent light source. The 
low energy electron beam can also produce the X-ray radiation. One of the methods for the 
low energy electron beam produced x-ray radiation is the Thomson scattering[2l]. When a laser 
beam ( with the frequency cjq )interacts the electron beam at an angle 0 , the upshift frequency 
of the laser as 


UJ = 


27^(1 — coscj)) 


1 _L ZQ. 
^ ^ 2 


-UJQ. 


( 1 ) 


where 7 is the lorentz factor, ao is the normalized laser strength. Using this mechanism, we 
can design the experiment scheme for x-ray radiation facility based on the 100 MeV electron 
beam from the laser plasma accelerator. Figure 22 shows the basic setup of the compact x-ray 
production facility. 
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Figure 22. the scheme for the x-ray experiment 


We can use the collimator to obtain the electron beam with the required energy spread. 
The typical parameters of the output electron bunch can be the charge of lOpC, the energy of 
lOOMeF, the energy spread of 5%, and the normalized emittance of O.lmmmrad respectively. 
We have used the laser beam of the wavelength of 3/im to do some laser electron beam scattering 
calculations. The output photon energy can be as high as 50KeV. 
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